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Abstract. We construct an exact relativistic cosmology in which an inhomogeneous 
but isotropic local region has fractal number counts and matches to a homogeneous 
background at a scale of the order of 10 2 Mpc. We show that Einstein's equations and 
the matching conditions imply either a nonlinear Hubble law or a very low large-scale 
CN ■ density. 
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1. Introduction 

O 

There is doubt and argument about whether the data on the galactic number count N 

5h . 

can support what is called a fractal structure, i.e., 



N<xy", (1) 

where < v < 3 is the fractal index and y is some distance measure (see |], ||, [| and 
references therein). However, there is evidence for such a pattern in the number counts 
out to a distance of the order of 10 2 /?, -1 Mpc |], @> El- 

Even this evidence is beset with problems because it is difficult to agree adequate 
statistics which are model independent p. For instance, Peebles (IJ, page 212) 
specifically excludes inhomogeneous spherically symmetric models in which we are at 
the centre. This may be reasonable but it serves to illustrate how aware one needs 
to be of the underlying model. It would help if good predictive dynamical models 
for the local distribution were available. The data could then be tested against these. 
Attempts in this direction are in their infancy but there are some that could yield a 
fractal number count [[J . At the same time, a number of workers are using recent large- 
scale structure observations and indirect evidence to argue that fractal-like clustering 
on small and possibly intermediate scales does not continue indefinitely, and that 
there is a transition to a near- homogeneous distribution on large scales H. This 
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is supported by the implications for large-scale structure of the near-isotropy of the 
cosmic microwave background radiation, since it has been shown that if all fundamental 
observers measure small anisotropy in the temperature, then the universe on large scales 
after last scattering is necessarily close to homogeneity J7|. 

Our study will consider the number counts as a function of the observer area distance 
(angular diameter distance) D or redshift z, down the past light cone of the observer. As 
a consequence, even in the spatially homogeneous Friedman-Lemaitre- Robertson- Walker 
(FLRW) case the number counts will be non-homogeneous. However the underlying 
spatial homogeneity of the FLRW geometry is reflected in the fact that the number 
count /redshift formula is precisely determined by the model ||. This formula specifically 
rules out a simple fractal distribution of the form flU). Note that our approach precludes 
direct comparison with work in which the fractal distribution is assumed to hold on 
spacelike surfaces (see for example ||). 

To analyse the problem theoretically we assume (a) that locally we are in a 
spherically symmetric region which can be modelled by a Lemaitre-Tolman-Bondi (LTB) 
geometry with a fractal number count of sources down our past light cone, and (b) that 
at some value of the observer area distance (of order lO 2 ^ 1 Mpc) the universe becomes 
homogeneous and can be modelled by an appropriate FLRW spacetime. We assume 
that the LTB region matches smoothly to the FLRW region, with no surface density 
layer or shell crossing [10, llj]. These assumptions are not unreasonable on the basis of 
current observations. Of course, we are only considering a very restricted form of fractal 
distribution, i.e. a single fractal in the radial direction which maintains the overall 
spherical symmetry. 

A similar model was investigated by Ribeiro ||12fl , who was the first to 
construct a relativistic fractal-count model. Ribeiro developed sophisticated numerical 
computations to analyze his models. He showed that parabolic and elliptic models gave 
a nonlinear Hubble law, while hyperbolic models were able to support a linear Hubble 
law. Our results are in agreement with these conclusions, but we believe that we have 
identified more clearly and simply the underlying reasons for these features. We have also 
pursued the implications of parabolic models that match to open FLRW models - i.e. 
that the large-scale density is extremely low. Ribeiro's work concentrated on numerical 
integration. Our approach is primarily analytical, and relies crucially on a systematic 
use of the central regularity conditions and the transition matching conditions. Thus 
our work may be seen as complementing Ribeiro's by bringing analytical insights to 
some of his results and extending aspects of his work. 

In principle, LTB models can accomodate a great variety of number count relations 
H O, 14|. However, the assumed number count relation may not be compatible with 



matching to the large-scale homogeneous universe, with other observations, or with 
regularity conditions on the observer's wordline. 
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The main result here is that if the fractal number counts are taken at face value and 
the universe can be modelled locally by an LTB geometry which matches to an FLRW 
background at some scale, then the universe either cannot have a linear Hubble law 
at low redshift (parabolic case), or it has a very low large-scale density (non-parabolic 
case). Perturbation of the model is unlikely to change this conclusion and increasing 
the scale of the fractal-count region (and there is some evidence for this) will exacerbate 
the low density problem. The low density problem identified by this model can only be 
avoided if the dark matter is strongly biased. 

It would be of interest to construct a linear perturbative model to confirm the above. 
The objective here, however, is to use an exactly isotropic non-perturbative model and 
to study the consequences flowing from the field equations and matching conditions. 



2. The cosmological model 



In this section equations are presented for the general LTB metric, which, together with 
a Kantowski-Sachs-type solution pTOfl , encompasses all regular, spherically symmetric 
dust models. These equations are used in section 3 to analyse the fractal number count 
subclass of LTB models. 

The metric must satisfy regularity conditions, including the central conditions as 
D — > 0. (See |J for the central conditions in general; a comprehensive analysis for 
LTB is given in [[KJ for central observers, and in |15[ for off-centre observers.) For this 
problem we use coordinates based on the observer's past light cones as described in 
H and usually called observational coordinates. Details of the constuction for the LTB 
model are given in il3fl . (See also where semi-isotropic observational generalizations 
of LTB models are investigated.) 

Explicitly the coordinates are {w,y,8,(f)}, where {w =constant} are the observer's 
past light cones, y is a distance from the observer along a past light ray and (6, 0) label 
the direction of observation. In these coordinates the metric of the LTB model is 

ds 2 = -A(w, y) 2 dw 2 + 2A(w, y)B(w, y)dwdy + C{w, yf (d6 2 + sin 2 0d<p 2 ) 

The Einstein field equations cannot be integrated explicitly in these coordinates 
13fl , but they can be in the more familiar (1+3) formalism fI7| . However the (1+3) 



coordinates are not directly observable. An advantage of the observational coordinates is 
that the null geodesic equations are already integrated, whereas they are not integrable 
in the (1+3) formalism. This explicit relation between the null geodesies and the 
coordinates facilitates the interpretations. We will use both formalisms below. The 
metric in (1+3) comoving coordinates has the form 

2 



ds 2 



-dt 2 + 



dR(r, t) 
dr 



dr 2 
1 - kf(r) 2 



+ R{r,t) 2 (d9 2 + sin 2 
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where / is arbitrary and relates to the 'total energy' of the LTB model. The cases 
k = 0, +1, — 1 correspond respectively to parabolic, elliptic or hyperbolic geometry of 
the {t = constant } hypersurfaces. If {w = wq} is the past light cone of observation, 
then 



D(y(z)) = C(w , y(z)) = R(r(z), T{r{z))) 



(2) 



where t = T(r) is the equation of a past light ray. 

In the observational coordinates the number count (total number of sources out to 
distance y) is given by || [13 



N{y) 



47r / n(wo, x)B(wo, x)D(x) dx . 
Jo 



where n(wo,y) is the number density of sources. The dust density is 



P 



mn . 



where m is the average galactic mass, which we assume to be constant. For simplicity, 
we omit evolution and selection effects. It would be possible to take some account of 
these effects when choosing n and m, but we have not done so. 

The total number count in the inhomogeneous region is given in terms of the Bondi 
mass function M(r) (which arises from integrating the field equations) by ||13 



N(y) 



m Jo 



M'(x) 



1 - kf(x) 



-1/2 



dx . 



(3) 



where M(y) is short-hand for M(r(y)). 

It is well known that the LTB metric is fully determined by three arbitrary functions: 
one can be removed by a coordinate choice but the other two have physical significance. 
In quadratures are given which determine the LTB functions directly from the 
observational data D(z) and N(z). This determination forms the corner-stone of the 
current analysis, since the LTB functions are needed to determine the dynamics and 
(crucially) the matching to the FLRW model. If the LTB functions are identified as the 
Bondi mass M, the energy function / and the big-bang time (3 (see |p.8|), then we can 



use the results of [13] to give these functions in terms of the observables as follows: 




2D (z 



where 



1 



Q(z 

By a rearrangement of ([!]) 

M(z 



ni 



x)N'(x) 



D(x) 



D'(x) 



dx . 



D{x)Q{x) 
(l+x) 2 



dx . 



(4) 



(5) 



m 



1 - kf(x) 



1/2 



N'(x) dx. 
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and 



(3{z) =t - D\x) \l - kf(x) 



1/2 



2M(x) 



D(x) 



kf(x) 



1/2' 



dx 



' -\[2D{zfM{z) 



-lll/2 



+ { 



M(z) [sinr(z)-r(z)]|/(z) 



(if k = 0) 

(if k = +l) 



(6) 



k M(z) [T(z) -sinhr(z)] |/(z)~ 3 | (if k = -1) 



where 



2 arcsin 



\D{z)f{zfM{z\ 



-i 



1/2 



r(z) 



2 arcsinh 



i J D(z)/(z) 2 M(z)- 1 



1/2 



(if k = +l) 



(if fc = -1] 



(7) 



An arbitrary constant of integration appearing in (6) has been identified as to, the time 
of observation (to is arbitrary since only the time elapsed after the big bang, t — /3, 
has physical significance). In these integrations the freedom in y on the light cone of 
observation has been used up by setting A(wo,y) = B(wo,y); this choice simplifies the 
calculation fl~3|l . 

Now it is evident that in the k = case only one of the functions N(z) and D(z) 
is arbitrary; for example N(z) may be obtained from (|]) by setting k = 0, once D(z) is 
known. This covariant constraint on the observational data may be found explicitly as 
follows. By © 



l + z 



DdQ fdN s 



-i 



and then 



mdD \dD 
with k = gives 



(8) 



Q— 

V dD 



-2m 



dN fd& 
dD \dD 



dD 



D 



1 + 



Q 2 



'l + z) 



Integrating by parts, and using the central conditions, we find that 



-2mN + 2m 



dN fdQ s 
dD [dD 



-i 



D 



m 2 Q 2 dN 
D [dD 



'dcr 

dD 



-2 



0. 



Solving this as a quadratic (and again imposing the central conditions to eliminate a 
spurious root), gives 

-i 



l_dQ 
QdD 



m dN 
DdD 



l2mN 
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Together with (|8|), this gives the solution 

. m dN 

1 + z — 1 — \ — — exp 




DdD 




(9) 



which appears to be a new result, and which is central to showing that the Hubble law 
is nonlinear (see below). Thus when k = 0, if N is known in terms of D, then (^]) gives 
z in terms of D. 

3. Fractal number counts 

The model we construct is isotropic about the observer and inhomogeneous out to a 
distance D = Dh of the order of 10 2 Mpc, corresponding to a redshift z below 10 _1 . 
The LTB spacetime metric in this local region satisfies the Einstein equations with a 
dust source which yields a power-law number count with fractal index, of the form ([!]). 
To ensure consistency with observational data, we also require a linear distance/redshift 
relation out to the homogeneity scale Dh, i.e. a linear Hubble law for small redshift. 
The Hubble constant is H = 100/i km s _1 Mpc -1 , with < h < 1. For the large-scale 
universe (D > Dh) , we require a homogeneous FLRW geometry. The Darmois matching 
conditions at Dh uniquely determine the parameters of this FLRW solution. 



In all cases the central conditions il3| demand that for small D 



N(D) = (j^& + 0{D*), (10) 

so that iV ~ D 3 as D — > 0. This means that the number count cannot be fractal 
(i.e. with v ^ 3) for D near zero. There is some minimum distance Df for fractal 
counts, below which N ~ D 3 . (Note that Df should be above the averaging scale which 
is implicit in a continuum dust model of galactic matter.) A simple model of fractal 
number counts that incorporates the limiting behaviour ( [TO] ) is the continuous power-law 
ansatz 

(^f) D 3 for D < D f 

N(D) = l ^ " ' (11) 



^) D ){w f T fOT Df<D<D h . 



(Matching conditions require iV to be at least continuous at D = Df, otherwise there 
would be a surface layer.) This model represents an 'instantaneous' transition from non- 
fractal to fractal counts at the minimum fractal distance Df, and N is not differentiable 
at Df. An alternative model, in which N is a smooth function of D, is 



which satisfies iV ~ D 3 for D < D f and N ~ D u ioi D ~> D f . 

Since has to be small, the region where D < Df can be treated as homogeneous 
and so Zf = z(Df) « H Df. At the homogeneity scale -D^, the number counts must 
match continuously to the FLRW number count relation N(D): 



N(D h ) = N{D h ) . 
The functions N{D) are known. For example, the k 



( 



mH - 



V 32vr 



1/3 



FLRW spacetime has |J 

2 



V 32 



7T 



(13) 



(14) 



5. 1. Parabolic fractal- count models 

Consider now the possibility of modelling the fractal-count region Df < D < Dh by a 
parabolic LTB solution (k = 0). Assuming that the core region < D < Df is also 
parabolic, Einstein's equation determines the Hubble constant in terms of the central 
density (as in a FLRW model) [|H 



r7rp 



We assume that the fractal number counts are modelled by the power-law ansatz ([TT1) . 
Then the matching of N at the homogeneity scale Dh, given by (|13D and fll4|) , implies 



u/3 



1 6 

- 8 









3-v 

^1 =o. 



(15) 



Thus the four parameters H , Df, D h and v are subject to the constraint (|15|) by virtue 
of number count continuity. If observations are used to determine D^ and v, then 
this constraint fixes the minimum fractal scale Df. 

The problem with the parabolic fractal-count models arises from the nonlinear 
behaviour of the redshift/ area distance relation for small D. From equations (Q) and 
([J), we find that (0) implies 

■D \ (3-i>)/2 



dD 



uHoDf ( 



(16) 



for small D. It follows that 

_ D (v-3)/2 



dz_ 
dD 



-D 



v-2 



for < v < 1 



for 1 < v < 3 



Thus after the initial linear behaviour up to Df (by construction), the redshift/distance 
graph begins immediately to curve downwards. This nonlinearity contradicts the well- 
established linear Hubble law on scales up to the order of 10 2 Mpc, and means that 
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parabolic fractal-count models are ruled out. Although we have used the power- law 



fractal count ansatz (11) to deduce this nonlinearity, it is clear that the feature will 
persist for any model that incorporates ([I]), since the argument depends only on the 
behaviour for small D. 

Finally, we note that exact expressions for z(D) may be obtained for any rational 
fractal index v, since in this case the quadrature in (0) may be performed exactly for 



the power- law relation flnp. 


For example, with v = 


-- f we 


find 


1+z = (1+Zf) 


'l-H D f (D/D f )^' 


2 

exp 


3H D f | 


1 - H D f 




(17) 



for Df < D < D] L . Nonlinearity for small D is readily confirmed by plotting the graph 

of m 



3.2. Non-parabolic fractal- count models 

The interpretation of observations cannot be realised if the fractal-count region is 
parabolic, since then the redshift/distance relation is nonlinear. For k — ±1 there 
is no such problem, since (as discussed in section 2) the observable functions N(z) and 
D(z) are independent. This allows us to circumvent the nonlinearity problem - but a 
new problem arises, i.e. the problem of very low density in the large-scale FLRW region 
arising from fractal number counts and the matching conditions. Intuitively, fractal 
number counts imply an under-density (since v < 3), so that k = — 1. Then matching 
conditions rule out k = +1 on large scales, and imply that the large-scale k = —1 FLRW 
model must have even lower density than the local LTB region. 

From (HI) the relationship between M and N is more complicated in LTB models 
with k 7^ 0. If k = +1, one obtains M/m < N, and conversely for k = —1. Note that 
regularity in the current context dictates that k must not increase with distance from 



the observer |T0[. Therefore we cannot match a k = — 1 LTB model to a k = +1 LTB 
or FLRW exterior. Furthermore, the matching conditions require continuity of M and 
kf\ 

Numerical integrations are simpler with differentiable functions. For this purpose 
we replace ( |TTD with the alternative smooth-transition ansatz ([12]) . We now focus on the 
dynamics and the implications of such a number count profile (with the other observa- 
tions) for the large-scale density. Using the formulas (f|)-(|7|) we have integrated the field 
equations with number count formula (|T2"D out to z = 0.07, where we assume the metric 
matches to the FLRW background. We take Df = 10 Mpc. We label the central (local) 
density parameter by fl c and the large-scale (background) density parameter by Qo. We 
take D(z) = Hq X z, which is the well established Hubble law for these redshifts. In the 
following tables the remaining parameters are varied in turn and the consequences for 
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the background density are shown. 



Solutions with h = 0.65, f2 c = 0.2 



fractal 
index v 


large-scale 
density Q 


1.0 


0.0002 


1.5 


0.001 


2.0 


0.008 


2.5 


0.04 


3.0 


0.2 



Solutions with h = 0.65, v = 1.5 



local 

density Q c 


large-scale 
density Q 


0.20 


0.001 


0.35 


0.0025 


0.50 


0.004 


0.75 


0.006 


1.00 


0.008 



Solutions with v = 1.5, f2 c = 0.2 



Hubble 

rate h 


large-scale 
density Q 


0.20 


0.00025 


0.40 


0.0007 


0.65 


0.001 


0.80 


0.002 


1.00 


0.0025 
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In all of these cases it was found that k had to be equal —1, corresponding to 
hyperbolic space geometry as one would expect with low densities. The background 
density is that of the uniquely defined FLRW model that matches at z = 0.07 (D = Dh) 
to the interior fractal-count LTB model (for more details of the matching calculation, 
see [TO]]). In all cases the significant fact is that the background density VLq is extremely 
low. Qualitatively, what happens is that the under-density implied by N ~ D u with 
v < 3 in the LTB region, is worsened by the fact that in the corresponding FLRW case, 
TV has a steep non-power- law gradient beyond 2 « 0.4 [J12| ]. 

It is of interest that the above numerical integrations confirmed that the no-shell- 



crossing conditions [11] were satisfied. The effect of the big-bang function (3 on the value 



of VLq is negligible at these redshifts. 
4. Conclusion 

The nonlinear Hubble law and low density problems of fractal-count universes, identified 
in this paper, apply to all the regular spherically symmetric dust spacetimes [[UJ . These 



are the spacetimes that can be constructed by piecing together regions with LTB metrics 
(including the homogeneous FLRW case), in which the matching satisfies the Darmois 
conditions and there are no surface layers or shell crossings. 

The nonlinear Hubble law at very low redshift rules out the parabolic models. 
However, ways to avoid the low density problem of non-parabolic models are still 
conceivable. There could be a surface layer or a cosmological constant, neither of which 
has been included in this analysis, and their effects are not fully known. The most 
compelling explanation, however, is the possibility of a significant bias in the data - the 
number counts of luminous matter may not trace the actual distribution of density. It 
could be the result of selection effects, evolution or the presence of dark matter. This 
provides an interesting slant on our result because it gives a different emphasis to the 
search for evidence. 
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